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FOREWORD 


The  science  of  underwater  warfare  has  progressed  rapidly  in 
recent  years,  particularly  in  the  development  of  torpedoes  of 
greatly  increased  range  and  speed,  and  homing  torpedoes  which  are 
required  to  run  complicated  trajectory  patterns.  This  advance  has 
been  accompanied  by  the  introduction  of  problems  in  the  control  of 
these  new  weapons.  Intuitive,  cut-and-try  methods  of  design,  which 
formerly  were  adequate  for  the  control  of  the  torpedo,  are  now  too 
expensive  and  time-consuming.  The  test  launching  of  a full-scale 
torpedo  is  an  operation  in  which  there  is  risk  of  damage,  or  loss 
of  the  missile.  Since  a torpedo  is  a costly  weapon  it  is  highly 
desirable  that  the  behavior  of  the  weapon  be  accurately  predicted 
before  it  is  ever  launched.  To  this  end  it  is  necessary  that  the 
equations  of  motion  and  the  mathematical  expression  of  the  laws  of 
motion  be  well  understood  and  be  expressed  in  usable  form. 

The  work  on  this  report  was  carried  on  under  Bureau  of  Ord- 
nance Task  Assignment  NOTS-C -6-257-16-54 . The  report  was  reviewed 
for  technical  adequacy  by  Milton  Plesset  of  the  California  In- 
stitute of  Technology  and  G.  V.  Schliestett  of  the  Naval  Ordnance 
Test  Station. 


N.  A.  RENZETTI,  Head 
Underwat  r Ordnance  Department 


Released  under 
the  authority  of: 

/ . ZCZ  ...  BRC'.'U 

Director 
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ABSTRACT 


Motion  equations  are  developed  for  a rigid  torpedo  of  con- 
stant mass  with  six  degrees  of  freedom.  It  is  assumed  that  the 
medium  is  at  rest  except  for  the  motion  caused  by  the  torpedo. 

The  mathematical  form  of  the  "mass  accession"  forces  is  derived 
from  potential  theory.  The  motion  equations  are  referred  to  body 
coordinates  In  their  development,  and  transformations  are  made  to 
inertial  coordinates.  Some  sources  of  hydrodynamic  coefficients 
are  discussed,  and  an  outline  of  the  methods  for  obtaining  them 
from  model  tests  is  presented.  Solutions  of  steady-state  equa- 
tions are  given,  as  well  as  a brief  explanation  of  the  analog 
computer  method  of  solving  the  motion  equations. 
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INTRODUCTION 


A great  deal  of  research  has  been  done  in  the  past  few  years 
In  the  analysis  of  the  motion  of  a torpedo.  The  problem  Of  un- 
controlled ship  and  torpedo  motions  was  investigated  by  Schlff  and 
Davidson  (Ref.  l)  and  by  Minorsky  (Ref.  2 and  3).  More  recently 
Bednarz  and  Harger  (Ref.  4)  studied  the  problem  with  the  purpose 
of  giving  a physical  insight  into  the  effects  of  the  coefficients 
of  the  simplified  motion  equations.  Triaxial  motion  equations 
were  developed  by  Pierce  and  Sepmeyer  (Ref.  5)  for  use  in- the 
Hydrodynamic  Simulator.  The  standardization  of  nomenclature  by 
the  Committee  on  Nomenclature  of  the  American  Towing  Tank  Confer- 
ence (1948)  was  an  important  step  forward  in  the  treatment  of 
motion  of  a submerged  body. 

In  the  past  the  most  serious  obstacle  to  the  analysis  of  tor- 
pedo motion  has  been  the  absence  of  information  concerning  the 
hydrodynamic  forces  acting  on  a torpedo.  To  eliminate  this  defi- 
ciency towing  tanks  and  water  tunnels  were  constructed.  Data 
obtained  in  model  tests  at  these  tunnels  yield  the  dimensionless 
coefficients  that  characterize  the  hydrodynamic  behavior  of  a 
torpedo . 

Techniques  and  facilities  for  the  solution  of  the  motion 
equations  have  been  expanded,  and  it  is  now  possible  to  solve  com- 
plex control  problems  with  relative  ease.  An  Important  facility 
used  in  the  analysis  of  torpedo  motion  is  the  electronic  analog 
computer.  With  the  aid  of  the  computer  linear  or  nonlinear  equa- 
tions are  solved  rapidly  and  accurately. 

In  this  report  the  motion  equations  are  developed  on  as  firm 
a theoretical  basis  as  present  knowledge  permits.  An  explanation 
of  the  methods  of  analysis  may  be  found  in  Ref.  6.  A discussion 
of  the  methods  by  which  hydrodynamic  coefficients  are  measured  in 
model  tests,  and  a very  brief  outline  of  the  analog  computer 
method  of  solution  of  the  motion  equations  are  included.  Since 
an  understanding  of  the  theory  of  torpedo  motion  must  precede 
applications  to  the  design  of  new  weapons,  this  report  Is  presented 
with  the  hope  that  investigators  in  the  fields  of  hydrodynamics  and 
torpedo  control  will  be  aided  in  understanding  the  present  state  of 
the  technique  and  be  stimulated  towards  its  advancement. 
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MOTION  EQUATIONS 


The  reaction  of  a torpedo  to  external  forces  Is  expressed  by 
the  fundamental  laws  of  dynamics.  In  the  development  of  the  motion 
equations  in  this  report,  it  has  been  necessary  to  make  assumptions 
about  the  nature  of  the  torpedo.  One  assumption  has  been  that  it 
is  a rigid  body.  The  torpedo  is  elastic  to  some  degree,  of  course, 
and  has  moving  parts —proper ties  which  will  probably  be  of  interest 
-n  future  studies.  At  present,  however,  it  is  felt  that  the 
assumption  of  rigidity  is  valid  for  the  applications  that  the  equa- 
tions of  this  report  have  in  view.  It  is  assumed,  moreover,  that 
the  torpedo  is  of  constant  mass.  This  is  a better  assumption  for 
electric  torpedoes  than  for  turbine  or  engine  driven  torpedoes, 
since  a considerable  quantity  of  fuel  is  consumed  in  the  latter. 
Usually  the  rate  of  fuel  consumption  is  so  slow  that  it  has  little 
effeot  on  the  trajectory.  In  a particular  problem,  however,  an 
investigation  should  be  made  as  to  the  length  of  the  trajectory 
for  which  changes  in  the  magnitude  of  the  inertia  of  the  torpedo 
may  be  neglected.  In  particular  problems,  moreover,  it  may  be 
necessary  to  modify  the  equations  as  they  are  given  here.  For 
example,  it  has  been  assumed  that  the  thrust  of  the  propulsion 
system  acts  along  the  longitudinal  axis  of  the  torpedo  without  re- 
sultant torque.  In  some  cases  it  may  be  necessary  to  add  addi- 
tional terms  to  the  equations  If  the  thrust  is  misaligned  or  if  an 
unbalanced  torque  is  present.  It  is  assumed  the  torpedo  is  fully 
wetted.  If  it  is  in  a cavltatlng  state  the  equations  given  here 
remain  applicable,  but  a modification  of  the  hydrodynamic  forces 
is  necessary  because  they  are  then  functionally  related  to  the 
cavitation  parameter. 


BASIC  VECTORIAL  EQUATIONS  OF  MOTION 

The  laws  of  motion  are  applied  to  a torpedo  under  the  assump- 
tion that  it  is  a rigid  body.  The  basic  equations  are 


(1) 

dQj) 

F «— — 

dt 

(2) 

dHjj 

dt 

where  F is  the  resultant  external  force  applied  to  the  torpedo 
body,  Gjj  is  its  linear  momentum,  L is  the  resultant  moment  acting 
on  the  torpedo,  and  Hj,  is  its  total  angular  momentum.  The  time 
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rate  of  change  of  the  torpedo  momenta  Gj>  and  must  be  with  re- 
spect to  an  inertial  reference  frame  for  the  application  of  the 
dynamic  equations.  A right-handed  rectangular  reference  frame 
( x0 , Hot  z0)  fixed  with  respect  to  the  earth  will  be  used  for  this 
purpose.  The  x0-y0  plane  is  tangent  to  the  surface  of  the  earth, 
and  the  z0  axis  is  vertically  downward.  Use  of  this  reference 
frame  as  an  Inertial  reference  frame  implies  the  assumption  that 
the  motion  of  the  earth  has  a negligible  effect  on  the  trajectory 
of  the  torpedo. 

An  underwater  missile  in  accelerated  motion  produces  accelera- 
tions in  the  flow  of  the  fluid  in  which  it  is  moving.  Consequently 
there  is  a transfer  of  kinetic  energy  to  the  fluid.  The  rate  of 
change  of  this  kinetic  energy,  and  therefore  the  force  producing 
it,  is  proportional  to  the  acceleration  of  the  missile.  Since  the 
inertial  reaction  of  the  missile  is  also  proportional  to  accelera- 
tion, the  missile  behaves  as  if  its  mass  were  increased.  This 
phenomenon  is  termed  “mass  accession".  It  will  be  assumed  in  this 
report  that  the  mathematical  form  of  the  "mass  accession"  forces 
is  given  by  the  theory  of  ideal  fluid  flow. 

Let  Fjl  and  Li  be  respectively  the  force  and  the  moment  (pre- 
dicted by  an  ideal  fluid)  on  the  torpedo,  and  let 

(3)  Pg  = P - P i 

00  Lg  * L - Lj 

The  force  and  moment  on  the  torpedo  predicted  from  potential  flow 
are  equal,  respectively,  to  the  negative  of  the  time  rates  of  change 
of  the  linear  momentum  and  the  angular  momentum  of  the  fluid  (see 
Appendix  B).  Thus 


(5) 

(6) 


d®f 


By  defining  the  system  of  body  and  ideal  fluid  momenta  as 


(7) 

Q = Of  + Gjj 

(8) 

H = Hj  + Bfa 
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the  basic  motion  equations  may  be  written 


dG 

(9) 

1 | 
n 

*°i 

dt 

dH 

(10) 

Lp  

dt 

It  is  advantageous  to  employ  a frame  of  reference  fixed  with 
respect  to  the  torpedo  body.  The  reference  frame  moves  with  the 
velocity  V of  the  torpedo  and  rotates  with  the  angular  velocity  co . 
In  this  reference  system  (see  Appendix  C) 

(11)  Fg  = G + cox  Q 

(12)  L2  = H + dOXH  + VXG 

o • 

where  G and  H represent  the  time  rates  of  change  of  G and  H as  seen 
from  the  moving  system. 

The  momenta  of  the  system  may  be  evaluated  from  the  total 
kinetic  energy  T of  the  system.  Let  T be  expressed  as  a function 
of  the  components  of  the  velocities  V and  in  the  body  reference 
frame.  Then,  letting 

(13)  V * iU  + Jy  + kw 
and 


(l4)  0)  = ip  + + kr 

the  momenta  are  given  by 


bT 

&T  *>T 

(15) 

G = i 

+ ~ — + k 

” dU 

„ dv  dw 

aT 

dT  bT 

(16) 

H = 1 

+ + k 

~ ip 

iq  br 

(see  Appendix  D) . The  total  kinetic  energy  T is  the  sum  of  the 
kinetic  energy  Tb  of  the  fluid  and  the  kinetic  energy  Tf  of  the 
torpedo  body.  These  are  derived  In  the  following  section. 
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KINETIC  ENERGY  OP  SYSTEM 


Kinetic  Energy  of  the  Torpedo  Body 

The  torpedo  body  will  be  considered  a rigid  aggregate  of  mass 
particles.  Let  mi  be  a representative  particle,  and  let  r*  be  the 
radius  vector  from  the  origin  of  body  coordinates  to  mi..  The 
velocity  of  m*  is 


(17)  s v + w x 

The  kinetic  energy  of  the  body  is  equal  to  the  sum  of  the 
kinetic  energies  of  the  individual  particles  of  mass,  and  is  given 
by 

(18)  Tt  . (l/2)Jmi|vi|2 

Let  r«  have  the  body  coordinates  Xj,  y«,  z*.  The  expansion  of 
Eq.  18  yields: 


(19)  2TJ,  * J mj/ljj2  + v2  + w2  + (yj2  + z^Jp2  + (x^  + z^Jq2 

+ (xi2  + yi2)  r2  + 2Uqzi  - 2Uryi  - 2qryizi  + 2vrx* 
- 2vpzi  - 2rpxizi  + 2wpyi  - 2wqxi  - 2pqxiyi| 

The  following  quantities  are  defined: 

(20) 


J mi  « m,  torpedo  mass 

CO 

•ri 

CO 

■P 

Jmi(yi2  + zi2)  = Ix 

fH 

Q> 

c 

J miyiZi  a 

lyz 

1 

^mjxi  * mxQ 

2»i(xi2  + zi2)  = Iy 

•n 

<M 

O 

Jmixizi  = 

Ixz 

<1h 

O 

^miyi  = myQ 

Jmi(xi2  + yi2)  = 1 2 

• 

} mixiyi  a 

Ixy 

•O 

O 

tg 

Jm±zi  = mzG 

& 


Equation  19  then  becomes 

(21)  2Tb  = mU2  + mv2  + mw2  + Ixp2  + Iyq2  + Izi*2 

+ 2mzgUq  - 2ray0Ur  - 2Iyzqr  + 2mxgvr 

- 2raZ(jvp  - 2Ixzrp  + 2myQwp  - 2cxQwq  - 2Ixypq 

The  origin  of  the  body  coordinate  system  is  usually  placed  on  the 
longitudinal  axis  of  the  torpedo  above  the  center  of  gravity.  The 
positive  x-axls  is  in  the  forward  direction  of  the  longitudinal 
axis,  and  the  positive  z-axis  is  vertically  downward  through  the 
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center  of  gravity.  The  choice  of  this  reference  frame  results  in 
the  simplification  of  the  expression  for  T^,  since  in  this  case 

*0  = ^0  = ^ = Iyz  = 0 
Kinetic  Energy  of  Ideal  Fluid 

Suppose  the  torpedo  to  be  moving  in  the  direction  of  its 
longitudinal  axis  with  the  velocity  U = U(t).  It  will  be  assumed 
that  the  flow  produced  is  irrotational  and  that  the  fluid  is  non- 
viscous  and  incompressible.  Let  the  velocity  field  of  the  flow  be 
q(x,y,z,t).  Because  the  flow  is  irrotational  a potential  function 
0(x,y,z,t)  exists  such  that 

(22)  a - - V0 

Since  the  fluid  is  incompressible,  the  divergence  of  the  velocity 
vanishes;  and  consequently 

(23)  dlv  £ - - V20  = 0 

There  can  be  no  flow  across  the  surface  of  the  torpedo.  Hence  the 
normal  component  of  a point  on  the  surface  of  the  torpedo  must 
equal  the  normal  component  of  the  fluid  velocity  at  that  point. 

Let  the  unit  normal  to  the  surface,  drawn  toward  the  fluid,  be 

(2t)  n = 1A  + + kV 

On  the  surface,  then, 

^0 

(25)  » AU 

bn 

It  is  assumed  that  the  flow  is  started  from  rest,  and  it  is  impos- 
sible that  finite  forces  acting  for  a finite  time  produce  a flow 
with  infinite  kinetic  energy.  Therefore,  the  velocity  of  the  flow 
must  vanish  at  an  infinite  distance  from  the  torpedo,  since  a 
finite  velocity  at  an  infinite  distance  would  imply  an  infinite 
kinetic  energy.  Hence  a potential  function  0 is  required  which 
satisfies  Eq.  25  on  the  surface  of  the  torpedo  and  whose  gradient 
vanishes  at  Infinity.  A solution  Is  sought  having  the  form 

(26)  0 = U(t)0x  (x,y,z) 

Since  ^01  = 0 with  ^0i/bn  = -X  on  the  surface  of  the  torpedo  and 
V0^  s 0 at  infinity,  the  function  0i  is  uniquely  determined.  Hence 
Eq.  26  is  the  solution  to  the  flow  problem,  since  an  irrotational 
flow  with  vanishing  divergence  is  uniquely  determined  by  its 
boundary  conditions. 
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Suppose  now  that  the  torpedo  is  rotating  about  Its  longitu- 
dinal axis  with  an  angular  velocity  p = p(t).  As  before,  a 
velocity  potential  exists  satisfying  Eq.  23-  Moreover,  at  a point 
(x,  y,  z)  of  the  surface  of  the  torpedo 

(27)  - — = P(t^  -Az) 

dn 

A solution  is  sought  of  the  form 

(28)  <p  = p(t y,  z) 

^ $4 

Since  v2^  = 0, — a Vy  -yttz  on  the  surface  of  the  torpedo,  and 

dn 

^ 4 * 0 at  infinity,  $4  is  uniquely  determined.  Hence  Eq.  28  is 
the  solution  to  this  flow  problem.  Several  velocity  distributions 
£l,  £2*  £3*  • • • may  added  together  to  obtain  another  velocity 
distribution.  The  velocity  distribution 

£ **  + £2  + £3  • • • 

is  said  to  be  the  superposition,  or  rather  the  result  of  the  super- 
position, of  the  velocity  distributions  £j,  £3*  • • • • It  is 

evident  that  if  the  divergence  or  the  rotation  of  each  of  the 
velocity  distributions  vanishes,  the  divergence  or  rotation  of 
their  superposition  vanishes  also. 

How  consider  the  flows  produced  by  motion  of  the  torpedo  in 
each  of  its  remaining  degrees  of  freedom.  For  motion  in  the  direc- 
tion of  the  y-axis  a potential  v'gfg,  and  for  motion  in  the  direction 
of  the  z-axls  a potential  wfo  is  obtained.  Rotary  motion  about  the 
y-axls  yields  a potential  q0j,  and  motion  about  the  z-axis  yields 
r0g.  The  potential  function  for  the  flow  produced  by  motion  of  the 
torpedo  in  its  six  degrees  of  freedom  is  obtained,  using  the  prin- 
ciple of  superposition,  as 

(29)  0 = XJ0i  + v02  + w03  + p04  + q05  + r0g 
Let  Tf  be  the  kinetic  energy  of  the  fluid.  Then 

(30)  Tf  « (1/2)^ (\70)2  dr 

where  the  integration  is  over  the  entire  volume  r of  the  fluid. 

The  Integral  may  be  transformed  by  Green's  theorem  to 
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(31) 


Tf  = -(l/2)/)/s^dS  -(l/2)/oLj2(720  dT 
ion  • 

r ^ of  _ 

= -(1/2)0 Js0*- — dS  since  = 0 
ion  • 


where  the  integration  is  over  the  torpedo  surfaces.  Substituting 
Eq.  29  into  Eq.  31  yields  a quadratic  form  in  the  torpedo  velocity 
components. 


(32)  2Tf 


anU2  + a12Uv  + ai3Uw  + 
+ a2]VU  + a 22^  + a^vw 

+ ag^wU  + a^wv  + a^W2 
+ ai^pU  + a42PV  + a^pw 
+ a51qTJ  + a52qv  + a^qw 
+ ag^rU  + ag^v  + ag^rw 


SifyUp  + a^sjUq  + 

+ a^vp  + a25vq 
+ a-^wp  + a^wq 
+ ai,i,p2  + ai^pq 
+ a5^qp  + a55q2 
+ agi^rp  + ag^rq 


ai6Ur 
+ a2gvr 

+ a^gwr 
+ a^gpr 
+ a5gqr 

+ a66r2 


where 

(33) 


aU  * Ws^i—  dS 
d n 


It  will  be  noted  that  Green’s  theorem  gives 


(3*0  aij  » aji 

Suppose  the  surface  of  the  torpedo  to  be  symmetric  with  respect  to 
the  x-z  plane  and  with  respect  to  the  x-y  plane.  Por  translatory 
motion  in  the  x-y  plane, 

(35)  2Tf  * anU2  + 2ai2Uv  + a22V2 


Because  of  the  symmetry  of  the  torpedo  the  kinetic  energy  must  be 
unchanged  if  v is  replaced  by  -v.  Hence  a12  a 0.  It  may  be 
similarly  shown  that  all  the  coefficients  of  cross-product  terms 
vanish  except  a^,  ag2,  a^g,  and  a^.  Equation  32  then  reduces  to 
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(36)  2Tf  = a^U2  + a2 2v2  + a^w2  + a^p2  + a^q2 

+ aggr2  + aajjwq  + 2a2gvr 

The  linear  momentum  and  the  angular  momentum  of  the  fluid  may  be 
obtained  from  Tf  as  Is  shown  in  Appendix  D. 

REACTION  OP  TORPEDO  AND  IDEAL  FLUID  SYSTEM 

The  total  kinetic  energy  of  the  system  composed  of  torpedo 
and  ideal  fluid  is  T * Tb  + Tf 

(37)  2T  = (m  + a1;l)U2  + (m  + a22)v2  + (m  + a33)w2  + (Ix  + a^Jp2 

+ (Iy  + 855) q2  + (Iz  + a 66 )*2  + 2mzoUq  - 21xzrp  - 2mzfivp 
+ 2agijwq  + 2a2gvr 

Define 

(38)  1 + an  » ^ 

b + a22  s m + B33  3 mip 
Ix  + ay,.  « Jx 

Iy  + a55  * Jy 

I z + a66  a Jz 

Then  Eq.  37  becomes 

(39)  2T  s m^U2  + mipv2  + mipw2  + Jxp2  + Jyq2  + J^r2  + 2mzQUq 

- 2Ixzrp  - 2mz(jvp  + 2a35wq  + 2a26vr 
The  components  of  momenta  defined  in  Eq.  15  and  Eq.  1 6 are 
(UO)  Qx  " mjjU  + mzqq 

Gy  a m>pv  - juzqP  + a2gr 

Gx  a BlipW  4 a-j^q 
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(40)  Continued  Hx  ■ Jxp  - Ixzr  - mzsv 

Hy  a Jyq  + mz6U  + a3jw 
Hz  = Jzr  - IXZP  + 

Let  ?2  have  the  components  Xg,  Yg,  Zg,  and  Lg  the  components 
K2,  M2 , N2«  The  basic  motion  equations  (Eq.  11  and  12)  written  in 
terms  of  their  components  become 

(41)  X2  = Gx  + qGz  - rGy 

Yg  B Gy  + rGx  - pGz 

Zg  a Gz  + pGy  - qGx 

Kg  = Hx  + qHz  - rHy  + vGz  - wGy 

Mg  » Hy  + rHx  - pHz  + wGx  - UGZ 

Ng  b Hz  + pHy  — qHx  + UGy  ~ vGx 

Substitution  of  Eq.  40  into  Eq.  4l  gives 

(42)  X2  » m^U  + mzgq  + m^wq  - vr)  + a^q^  + mzGpr  - agfj**2 

• • • 

Y2  = mrj«v  - mzop  + aggr  + mLUr  + mzgqr  - m^wp  - a^pq 

Zg  b mTw  + a^q  + mijvp  - mzg(p2  + q2)  + aggpr  - m^Uq 
• • • 

K2  = JXP  - Ixzr  - mzGv  + (Jz  - Jy)qr  - Ixzpq  + (a26  + a35)vq 

- mzGUr  - (a35  + a2g)  wr  + mzQwp 

• • « Q Q 

Mg  a Jyq  + mz6U  + a^w  + (Jx  - Jz)pr  + Ixz(p  - r ) - mzQvr 

- aggvp  + (mL  - m^JU  + mzQwq  - a^Uq 

Hg  = Jzr  - Ixzp  + aggv  + (Jy  - Jx)pq  + a35wp  + Ixzqr 
+ (mip  - mL)Uv  + aggUr 
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These,  then,  are  the  basic  motion  equations  for  a torpedo  referred 
to  a set  of  axes  fixed  with  respect  to  the  torpedo,  the  positive 
x-axis  in  the  forward  direction  of  the  longitudinal  axis,  and  the 
z-axis  vertically  downward  through  the  center  of  gravity  of  the 
torpedo.  Assumptions  under  which  they  have  been  derived  are 

1.  The  torpedo  is  a rigid  body  of  constant  mass,  symmetric 
with  respect  to  a plane  through  its  longitudinal  axis. 

2.  Motion  of  the  earth  has  negligible  effect  on  the  trajec- 
tory of  the  torpedo. 

3.  The  medium  is  infinite  in  extent  and  at  rest  except  for 
the  flow  produced  by  the  motion  of  the  torpedo. 

The  torpedo  is  fully  wetted. 

5.  Mass  accession  forces  are  formally  given  by  the  theory  of 
ideal  fluid  flow. 


EXTERNAL  FORCES  AND  MOMENTS  ACTING  ON  THE  TORPEDO 

The  external  forces  and  moments  acting  on  the  torpedo  are 
those  caused  by  gravity  and  the  propulsion  system,  and  those  pro- 
duced by  hydrodynamic  and  hydrostatic  pressures.  The  net  force 
produced  by  gravity  is  the  weight  of  the  torpedo  acting  vertically 
downward  at  the  center  of  gravity.  The  net  force  of  the  hydrostatic 
pressures  is  a buoyant  force  acting  at  the  center  of  buoyancy  of 
the  torpedo.  Resolution  of  these  forces  and  the  moments  produced 
by  them  onto  body  coordinates  is  given  in  Appendix  E.  The  thrust 
of  the  propulsion  system  will  be  assumed  to  act  along  the  longitu- 
dinal axis  of  the  torpedo  with  no  resultant  torque. 

Hydrodynamic  forces  predicted  by  the  theory  of  ideal  fluid 
were  discussed  above.  This  theory  predicts  no  lift  or  drag  on  the 
torpedo.  Since  a torpedo  does,  in  fact,  experience  lift  and  drag 
forces,  they  are  ascribed  to  deviation  of  the  fluid  motion  from 
potential  flow  because  of  the  viscosity  of  the  fluid.  At  very  high 
Reynolds  numbers  this  deviation  from  potential  flow  will  take  place 
in  a thin  layer  in  the  neglhborhood  of  the  surface  of  the  torpedo 
(see  Ref.  7)*  Particles  of  fluid  at  the  surface  of  the  torpedo 
adhere  firmly  to  it  so  that  not  only  the  normal  component  of  the 
fluid  velocity  at  the  surface  but  also  the  tangential  velocity  is 
equal  to  that  of  the  surface.  Outside  the  boundary  layer  a poten- 
tial flow  will  exist.  Let  it  be  assumed  that  this  potential  is  the 
same  as  Eq.  29* 
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The  Navler-Stokes  equations  for  viscous  fluid  flow  give  the 
force  per  unit  volume  acting  on  the  fluid  as 

(43)  f = - VP  + kV2^ 


where  P is  the  pressure,  £ is  the  velocity  field,  and  k is  the 
coefficient  of  viscosity  of  the  fluid.  Integrating  over  the  entire 
volume  of  the  fluid  gives  the  total  force  acting  on  the  fluid  as 


(44) 


Ij.1  dT  = -fTv p d r + kjrv2a  dr  = JgPn  ds 


where  the  surface  integrals  are  over  the  surface  of  the  torpedo, 
and  n Is  the  unit  normal  to  the  surface  projecting  "into  the  fluid. 
It  ilT  assumed  that  at  the  upper  limit  of  the  boundary  layer 
£ = -v0-  At  a point  of  the  surface  whose  radius  vector  from  the 
origin  of  body  coordinates  is  r,  the  velocity  is  V + w x r.  Let 
the  thickness  of  the  boundary  layer  be  denoted  bylf . Then, 
approximately. 


(45) 


— » -(V0  + + w X r)^"1 


The  difference  In  pressure  between  the  inner  and  outer  surfaces  of 
the  boundary  layer  Is  small  (Ref.  7)*  Hence,  the  additional  force 
on  the  torpedo  because  of  the  viscosity  of  the  f?uid  is  approxi- 
mately 

(46)  -/s(V0+  V +&J  x rjktf”*  dS 


Since 


(V  +w  X r)-n  « -(V0)»n  at  the  surface, 

-V0-  (V  + W x p)  represents  the  relative  velocity  at  which  the 
potential  flow  Ts  sliding  over  the  surface.  This  is  a linear 
function  of  the  torpedo  velocity  components.  If  cT  were  independent 
of  the  velocity  components,  then,  the  viscosity-induced  force  would 
be  al&o  a linear  function  of  the  torpedo-velocity  components.  This, 
however,  is  not  the  case.  The  boundary-layer  thickness  is  a func- 
tion of  Reynolds  number  and  the  form  of  the  surface.  Moreover, 
the  boundary  layer  may  become  quite  thick  and  separate  toward  the 
after  end  of  the  torpedo,  and  the  analysis  given  above  is  then  not 
applicable.  It  is  at  least  reasonable  to  assume,  however,  that 
the  viscosity  forces  and  moments  are  functions  only  of  the  torpedo 
velocity  components . 

The  remaining  hydrodynamic  forces  to  be  taken  into  considera- 
tion are  those  produced  by  deflections  of  the  control  surfaces. 
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These  forces  depend  primarily  on  the  magnitude  of  the  deflection 
of  the  control  surface  and  not  on  the  rate  In  control  systems  used 
at  present.  It  will  be  assumed,  then,  that  the  forces  and  moments 
produced  by  control  surface  deflections  are  functions  of  the 
magnitudes  of  the  deflections,  d"e  and  cfr. 

The  forces  and  moments  of  Eq.  11  and  12  may  therefore  be 
written 

(47)  Pg  = v,  w,  p,  q,  r,  cTe,  <fr)  + B + W + T 

Lg  = Lg(U,  v,  w,  p,  q,  r,  <fe,  <fr)  + rBxB+rfixW 

where  B is  the  buoyant  force,  W is  the  weight  of  the  torpedo,  T Is 
the  thrust,  rg  Is  the  radius  vector  to  the  center  of  buoyancy,  and 
tq  is  the  radius  vector  to  the  center  of  gravity. 

Let  the  components  of  Fg  be  Xg,  Yg,  Zg  and  the  components  of 

L be  Kg,  Mg,  Ng.  These  forces  and  moments  are  usually  determined 
In  model  studies  (see  section  entitled  "Sources  of  Hydrodynamic 
Coefficients").  Results  obtained  from  these  studies  show  that, 
for  most  torpedoes,  the  hydrodynamic  forces  and  moments  are 
approximately  linear  functions  of  angle  of  attack  and  turning  rate 
over  the  normal  operating  range.  Forces  and  moments  produced  by 
control-surface  deflections  are  also  linear  over  a wide  range  for 
most  torpedoes.  There  Is  some  Justification,  therefore,  for 
expanding  the  components  of  Fg  and  Lg  in  a Taylor  series  about 
U c U0,  v s w » p a q = r » rf'j  * = 0,  where  U0  Is  the  operating 

velocity  of  the  torpedo,  and  neglecting  all  but  first  order  terms. 
Hence,  letting  XgQ  be  Xg  evaluated  at  the  point  about  which  the 
series  expansion  Is  made,  and  similarly  with  the  other  components, 
and  letting  D » D0  + u, 

(48)  Xg  a XgQ  + XgyU  + XgyV  + Xgww  + Xgpp  + Xgqq 

+ Xgjj?  + Xg^cy;  + Xg^yj, 

Yg  3 Y20  + Y2uU  + Y2vV  + Y2»w  + Ygpp  + Y2qq 
+ YgpX*  ♦ Yg^/e  ♦ *2#* 

Z2  » Z20  + Z2uU  + z2vV  + Z2w*  + z2pP  + Z2qq 
Z2rr  + Z2cfe^e  + Z2^r 
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(48)  Contd . K2  - KgQ  + ^2u'1  x2w'*  x2q^ 

+ K^r  + K2(^<4  + K2(/r  4 
Mg  — Mgo  + MgyU  + MgyV  + MgyW  + Mgpp  + Mgqq 


+ Mgj.r  + ^2cf/e  + M2  cfifi 
N2  a ®20  + ^2uu  + ^2vv  + n2ww  + ^2pP  + N2q^ 


+ NgpT  + ^2cP/q  + *2dfi 


where  the  partial  derivatives  are  evaluated  at 

UsVsIfap.qsra^a^sO 

Because  of  symmetry  with  respect  to  the  x-z  plane;  and  because  of 
the  point  about  which  the  Taylor  expansion  is  made,  the  following 
partial  derivatives  and  component  values  at  the  point  of  expansion 
vanish: 


x2v  “ x2p  * x2r  3 x2q  ■ x2w  “ x2 tfe  3 x2/r  * ® 

Y2u  " Y2w  * y2q  3 y2  de  3 Y20  3 0 

Z2u  ■ Zgy  s Zgp  a Zgj,  a Zfcfr  a ZgQ  a 0 

*2u  * *2w  3 K2q  3 K2je  3 *20  3 0 

Mgu  ■ Mgy  8 Hgp  a Mgj,  a MgJ'p  a MgQ  a 0 

N2u  3 N2w  = H2q  3 N2<fe  3 N20  3 0 

Using  the  resolution  of  buoyancy  and  gravity  forces  and  mo- 
ments as  given  In  Appendix  E and  the  linearized  external  forces 
(Eq.  48)  the  equations  of  motion  (Eq.  42)  become 

(49)  X20  + *2uu  + 1 _ (w  “ B)  sin  ® * *LU  + mzGQ  + ®t(w9  - vr) 

+ a35q2  + mzQpr  - a26r2 
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(^9)  Y2vv  + Y2pP  + Y2rr  + ^/r^r  + (W  - B)  sin  0 cos  Q 
Contd. 

• * t 

= mipv  - mzGp  + a2gr  + m^Ur  + razQqr  - m^wp  - a3ejpq 
Z2ww  + Z2qq  + Z Q(fe^e  + (W  - B)  cos  0 cos  6 

= m^w  + a^q  + m^vp  - mzG(p2  + q2)  + a2gpr  - mLUq 

k2vv  + ^pP  + K2rr  + K2(/'p<:^*  " Wzg  sin  0 cos  ® 

- JXP  - Ixzr  - mzQ^  + (Jz  - Jy)  qi*  - Ixzpq  + (a2g  + a35)vq 

- mzgUr  - (a^j  + a2g)wr  + mzGwp 

m2ww  + M2qq  + ^2cT  <4  + BxB  cos  0 cos  ® " WzQ  sin  0 
© 

» Jyq  + mzGu  + a35w  + (Jx  - Jz)pr  + Ixz  (p2  - r2) 

- mz0vr  - a2gvp  + (mL  - aupjuw  + mzffwq  - a^Uq 
N2vv  + N2pp  + N^r  + ^2(tr^r  “ BxB  sln  0 cos  6 

“ Jzr  - IxzP  + a2gv  + (Jy  - Jx)pq  + a35wp  + Ixzqr 
+ (mip  - m^JUv  + ajjgUr 

The  motion  equations  (Eq.  49)  will  be  rewritten  in  terms  of 
angle  of  attack  a and  angle  of  sideslip  p.  By  definition 

(50)  sin  a = — 

V 

v 

sin  B = 

V 

It  will  be  assumed  that  a and  3 are  small  angles  so  that  the 
sine  of  the  angle  is  approximated  by  the  angle.  Moreover,  a change 
of  notation  will  be  Introduced  at  this  point.  Let 
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*o-*20 

Xu  = Xgu 

Ve  - *We 

vr  = k24 

Kp  = -KgvV 

Kp  «=  K2p 
Kr  = Kgjt 

Np  = -M2vV  + (tut  - miJV2 

Np  = »2p 

Nr  =»  Ngj.  - a2gV 

MP  = a26V 

%r  - N2<fp 


Yp  = -Y^V 

Yp  --  Ygp 

Yp  = Y^ 

Yr  " “a26 

vr  - Y2^r 

Ma  = M2wV  + <raT  - mL>v2 
Mq  * M2q  + a35v 
Ma  - ~a35v 
M<fe  = M2^e 

za  = z2wv 

zq  * z2q 

Zq  s "a35 

H = Z2< 


With  this  change  of  notation  the  equations  of  motion  (Eq„  49)  become 
(51)  XQ  + Xuu  + T - (W  - B)  sin  G 

= nijU  + mzQq  + mTV(aq  + pr)  + mzgpr  “Zqq2  + Vp?2 

Yp3  + YpP  + Yrr  + Y£r  + Y^  cfj,  + (W  - B cos  & sin  <P 
* • 

= - mjVp  - mzgp  + mLUr  + mzgqr  + Zqpq  - SbpVap 
Zaa  + Zqq  + Zqq  + <*e  + (w  ” B)  cos  0 cos  6 

- nkjVa  - nupVpp  - mzg(p2  + q2)  - Y£  pr  - mLUq 


NAVORD  REPORT  20S 

(51)  Kflp  + KqP  + K«r  + Ylf*  dL  - WzQ  sin  0 cos  9 
Contd.  p v r 

= JXP  - Ixz3^  + mzGv3  + (Jz  “ Jy)qr  " ^zPI 

+ V(Yp  + Zg)  (3q  + or)  + mzGUr  + raZQVap 

Maa  + Mqq  + M^a,  + M^c^  + Bxg  cos  0 cos  9 - Wzq  sin  9 

- Jyq  + mzGu  + (Jx  - Jz)pr  + Ixz(p2  - r2) 

+ mzGV3r  + Np3p  + mz0Vaq 

Np3  + NpP  + Npr  + N33  + Nc f^<PT  - BxB  sin  0 cos  9 
• • 

= Jzr  - XxzP  + (Jv  " ‘ *&» P + Ixz(lr 


TRANSFORMATIONS  OF  MOTION  EQUATIONS 

The  motion  equations  may  be  written  in  terms  of  the  inertial 
angular  position  of  the  torpedo, 0,  9,  <p,  (Ref.  4).  The  com- 
ponents of  the  rotational  velocity  of  the  torpedo,  p,  q,  and  r,  in 
body  coordinates  are  related  to 0,0,  (p  by 

(52)  p 3 (p  - "Y  sin  0 

m m 

q = ]^  cos  Q sin  0+0  cos  (p 
• • 

r = 0 cos  0 cos  <p  - Q sin  (p 

(see .Appendix  F).  Components  of  the  angular  acceleration,  p,  q, 
and  r,  are  subsequently  found  to  be 

(53)  P « (P  - Y sin  0 -00-cos  9 

q a 0cos  9 sin  (p  -00  sin  9 sin  (p  +(f(p  cos  <p  cos  6 
• • • * 

+ 6 cos  <p  -&p  sin  (p 

, (•  • • • » 
r = y cos  9 cos  (p  - y& sin  9 cos  <p  -Tf/tp  cos  6 sin  <p 

••  • . 

- & sin  (p  -dp  cos  <p 
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In  addition  to  the  assumptions  already  made  it.  will  be  assumed 
that  the  pitch  angle  $ is  small  so  that  approximately,  sin  9 ~ 
cos  0=1.  Using  Eq.  52  and  53  the  motion  equations  (Eq.  51)  are 
expressed  In  terms  of  Inertial  angles  as 

(54)  X0  + V + Tx  - (W  - B)0 

= m^u  + mzG(^sin  Cp  - ye&  sin  (f)  + ypoos  cp  + 9 cos  <p 

- Qp  sin  (p ) + mip V \y(y  sin  p + 9 cos  Ip) 

+ $(yr  cos  cp  - Q sin  - Z sin  cp  + 9 cos  p)2 
+ Yr(f  cos  (p  - 6 sin  Cp )2  + td.Zq((P  ~ Y&){Y  003  9) 

YpP  + Yp(^  - ye)  + Y r(Y  cos  cp  - 6 Sin  <p)  + Yp(j^  eos  cp 

- yed  cos  <p  - }/<p  sin  cp  - 9 sin  cp  - Qp  cos  p) 

+ Y (f  (Tt  + (W  - B)  sin  Cp  - -m«pYp  - mzQ(^  -ye  -Oy) 

+ mjU(^  cos  <p  - G sin  cp ) 

+ mzG(^  cos  cp  - 9 sin  <p)(y  sin  cp  + 9 cos  <p) 

- xs^a.{<}>  -jre)  + -ye)<y  sin  cp  + e ccs  <p) 

Zaa  + Zqiy  sin  tp  + 9 cos  (p)  + Z^y  sin  <p  - yB9  sir.  P 

+ yp  cos  cp  + 6 cos  cp  - £tp  sin  p)  + + (tf  ~ B)  e5s  9 

= mpVd  - nupVp  ( tp  -ye)  - mzG(<?  -fG)2 

- mzG($/  sin  <p  + 9 cos  p)2  - Y«,(^  - Y$)(p  cos  9? 

- 0 sin  0?)  - nuU(^  sin  p + 6 cos  9?) 

Kp3  + Kp(^  - ye)  + Kp(^  cos  cp  - 6 sin  p)  + Kj  - Ws&  3 in  p 

= jx$  - ^0  - 0$  - iX2(^  cos  9 ~ yee  cos  <p  ~yp  am  cp 
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(5*0  - & sin  <p  -9p  cos  $ + mz5v£ 

Contd. 

+ (Jz  - Jy)(^/cos  <p  - Q sin  $Kj£sin  <p  + 6 cos  <p) 

- Ixz6*>  sin  <p  + 0 cos  (p)  + V(Y*  + Z^)p 

(,^sin  cp  + 6 cos  00)  + 'mzGVa((^  - 'j/'d) 

+ V(Y^  + Z^JaC^cos  (p  -6  sin  <p) 

- mzgV^  cos  <p  - & sin  <p) 

• • 

Maa  + Mq(^sin  <p  + & cos  0?)  + M^a,  + M^c/^  + BxB  cos  cp  - Wz 

= Jy(^  sin  00  -^90  sin  cp-^yp  cos  <p  + 9 cos  - & sin  <p) 

+ ijizqU  + (Jx  - Jz){^  -^0)(^cos  <p  -Q  sin  cp) 

+ Ixz(<?  -$0)2  - IXZ(^C0S  0?  -0  sin  <p)2  + Np3(00  -ye) 

+ mzGVa(^sin  00  + <9  cos  0?)  + mzGVp(^  cos  0?  -<6  sin  <p) 
%p  + np (cp  -ye)  + Hr( y cos  <p  - b sin  cp)  + n^p 
+ N^dj.  - BxB  sin  <p 

= Jz(^cos  p -free  cos  & -ifop  3in  00-^  sin 
-&p  cos  p)  - htf  fa  -if) 

+ (Jy  - j x)0?  -'p‘9)\Y sin  <p  + 0 cos  00) 

- M^(0P  -y0)o-  + Ixz(^cos  (p  - 6 sin  0>)(/ sin  cp 
+ 5 cos  $ 

These  equations  are  greatly  simplified  if  terms  involving 
products  of  velocities  may  be  neglected.  Dropping  such  terms, 

Eq.  51*  becomes 
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(55)  x0  + XpU  + Tx  - (W  - B)0  = mjii  + mzg^sin  p + 9 cos  cfi 
YpP  + Yp(^  - f(9)  + Yr(^cos  p -Q  sin  p) 

+ Yp(^  cos  p - 9 sin  + Y^cfj,  + (W  - B)  sin  ^ 

= -nifpVfS  - m Zq$  -fyd)  + mLV(^'eos  g. 0-9  sin  p) 

Zaa  + Zg(^  sin  p + Q cos  <p)  + Zq($*  sin  (p  + 0 cos  p) 

+ + (w  “ B)  cos  p 

- mipVa  - mjV(^r  sin  p + 6 cos  p) 

KoP  + Kp(0>  - ye)  + YLpOfr  cos  p - 6 sin  $ + K,/cC,  - WsG  sin  ^ 
p v r 

= Jx(#  -$fe)  - Ixz(y  cos  0?  - e sin  p)  + mzGVp 
+ mzQV(y  cos  p - 6 sin  pf) 

Maa  + Mq(^  sin  p + 9 cos  00)  + M^a,  + + BxB  cos  p - WzG0 

= Jy(y  sin  p + 6 cos  p)  + razGu 
npp  + np (p  - ye)  + Nr(^  cos  p - e sin  p 
+ Njjp  + Kffyfp  - BxB  sin  p 

* Jz(^cos00  -6  sin  0 - Ixz('p-fi9) 
where  it  has  been  assumed  that  U0  S5£  V. 

If  <p  s Ixz  = Zq  = Yp  = Np  = 0,  and  the  forward  velocity  is 

constant,  Eq.  55  nay  be  reduced  to  two  sets  of  equations,  the  yaw 
equations  and  the  pitch  equations.  The  yaw  equations  are 

(56)  YoP  + Yrf  + l±f+  V cTr  = -JftpVp  + mLV^ 

H r 

NpP  + Nr?r  + NpP  + V/r  - hf 
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The  pitch  equations  are 

(57)  Zaa  + Zq0  + Zg$+  Z^eTe  + (W  - B)  - - mLV0 

Maa  + + Ve<^e  + BxB  = Jy® 


It  is  sometimes  convenient  to  change  the  origin  of  the  body- 
fixed  coordinate  system.  Suppose  the  center  of  gravity  is  to  be 
shifted,  for  example.  It  is  then  desirable  to  write  the  equations 
with  the  new  center  of  gravity  as  the  origin.  Hydrodynamic  forces, 
which  depend  only  on  the  exterior  shape  of  the  body,  are  unchanged. 
There  will  be  a change  in  the  hydrodynamic  coefficients,  however, 
because  of  the  change  in  angle  of  attack  at  the  new  origin.  Sup- 
pose the  origin  to  be  shifted  forward  on  the  longitudinal  axis  a 
distance  A,  and  consider  the  hydrodynamic  force  and  moment  to  be 

(58)  Y « Y2vv  + Y^r  + Y2^4 

» - n2,v  + or 


the  lateral  velocity  at  the  new  origin  is 

(59)  v*  = v + Ar 
The  moment  around  the  new  origin  is 

(60)  N*  = N - YA 


In  terms  of  v*,  then,  the  force  and  moment  at  the  new  origin  are 

(61)  Y . Y^v*  + (Yj*  - AY^Jr  + Y 

N*  = (Ngy  - AY^Jv*  + [Ngp  - AY^  - A(N2v  - AY^jJ  r 


Hence 

(62) 


+ (N2ef  - &SSK 


Ygv*  - ?2v 


Y * rn  Y-  - AY- 
2r  2r  2v 


W = *afr 
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(62)  N2v*  = N2v  - AY2v 
Contd . 

Ngj*  a Mgp  - AY^  - A(N2v  - AY2v) 

N2dr*  = N2cfr  - ^2</r 

Substitution  of  Eq.  59  into  the  expression  for  % (Eq.  36)  gives 

(63)  a gg*  = a66  - 2Aa2g  + a22A2 

a26*  " a26  - a22* 

The  yaw  equations  (Eq.  56)  become 

(6M  Yp*p*  + Yp*^+  Yp*^+  Y <fr*4  = -n^VS*  + mLV^ 

Np*p*  + Hp*y-+  Np*p*  + = J2*j/ 

where 

(65)  Yp*  » Yp 

V = 

y <rr*  « Y^r 

Yr*  * 'a26* 

Np*  = Np  - AYp 

Np*  = Nr  - AYr  + ~(Np  - AYp)  + anAV 

Np*  = a26*V 
N/p*  » N/p  - AY/P 

V = V + a66* 


Similar  transformations  may  be  carried  out  for  the  pitch  equations 
(Eq.  57) • The  new  pitch  equations  are 
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(66)  Z^*a*  + Zq*0  + Zq*0  + 2^*c£  + (V  - B)  > n^Va*  - mLV6 

Ma*a*  + MqS  + M£*a*  + M^*oQ  + BxB  = Jy*©‘ 


where 

(67) 


V - z, 


a 

Z * = Z_  + ^Z_ 
q q y a 


Z<C  " H 


zq*  = “a35* 


M * = M + AZ 
a a a 


Mq*  = Hq  + AZq  + |(Ma  + AZa)  + anVA 


Ma*  = -*35^ 


ty*  = + AZ^ 

e e e 


Jy*  53  Jy*  + ®55* 


a55*  * a55  + ^35  + a3J 


a35*  = a35  + a33* 


SOURCES  OF  HYDRODYNAMIC  COEFFICIENTS 


MODEL  TESTS  AS  SOURCES  OF  HYDRODYNAMIC  COEFFICIENTS 

Model  tests  may  be  classified  Into  three  types:  static  tests, 
rota ting -am  tests,  and  forced-oscillation  tests. 

Static  Tests 


In  static  tests  the  model  nay  be  towed  through  a tank  or 
placed  in  a water  tunnel  at  a fixed  angle  of  attack.  The  resultant 
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forces  and  moments  are  then  measured.  Suppose  the  model  to  be 
towed  through  a tank  at  a fixed  angle  of  sideslip  p,  and  let  the 
measured  force  and  moment  be,  respectively,  YM  and  NM.  Prom  Eq.56 

(68)  Ym  + Ypp  + V/r  = 0 

NM  + + V/r  = 0 

By  repeating  the  experiment  at  different  sideslip  angles,  rudder 
settings,  and  velocities  Yg,  Y^ , N«,  and  V may  be  determined  as 
functions  of  p,  efj>,  and  velocity.  By  rollingrthe  model  through 
90  degrees,  the  coefficients  Za,  Z , M®,  and  may  be  deter- 
mined in  a similar  manner.  e 

The  drag  Xg  may  be  determined  by  measuring  the  force  component 
along  the  longitudinal  axis.  It  is  found  that  for  the  range  of 
velocities  in  which  most  torpedoes  operate  X0  may  be  expressed  in 
terms  of  a dimensionless  coefficient  X^  and  the  velocity  as 

(69)  _ X0  = lMxJv2 

where  A is  a characteristic  area. 

The  coefficient  Xu  may  be  obtained  from  Eq.  2 by  differen- 
tiating It  with  respect  to  V.  Thus 

(70)  Xu  = /°AXqV 


Rota ting -Arm  Tests 

In  rotating -arm  tests  the  model  is  towed  In  a circular  path 
at  fixed  angular  velocities  and  at  a fixed  angle  of  attack.  For 
this  condition  the  motion  equations  (Eq.  56)  become 

(71)  Ym  + Ypp  + (Yr  - mLV)^+  Y<//r  = 0 

Nm  + Npp  + Nrr+  V/r  = 0 

Prom  these  equations  the  coefficients  Yr  and  Nr  may  be  deter- 
mined. The  coefficients  Y^  and  cannot  be  measured  with  this 
type  of  test  since  an  accelerated  motion  is  required.  Measurement 
of  these  coefficients  may  be  accomplished  by  forced  oscillation 
tests . 
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The  rotating-arm  test  has  the  advantage,  however,  of  making 
possible  the  determination  of  the  nonlinearity  of  the  hydrodynamic 
coefficients.  For  example,  M_  may  be  determined  as  a nonlinear 
function  of  a by  measurement  of  the  moment  M over  a range  of  angle 
of  attack  a.  For  a more  complete  discussion  of  the  nonlinearities 
of  the  coefficients  see  Ref.  4. 

Forced-Oscillation  Tests 


In  one  such  test  the  model  is  supported  by  means  of  a shaft 
In  a water  tunnel.  The  shaft  Is  made  to  oscillate  through  the 
application  of  a sinusoidal  torque  applied  through  a spring. 
Measurements  of  amplitude  and  phase  then  permit  the  computation 
of  the  hydrodynamic  coefficients. 

The  model  Is  placed  in  the  water  tunnel  as  3hown  In  Fig.  1. 


Spring 


FIG.  1 

Let 

« spring  constant  of  main  drive  spring 
K2  = spring  constant  of  shaft 
y1  b input  displacement  angle 
>2  8 output  displacement  angle 

Define  Xi  and  positive  In  the  same  sense  as  the  yaw  angle  y . 

Two  or  three  support  positions  may  be  used  In  the  test.  The  hydro- 
dynamic  coefflclehts  will  be  defined  about  one  of  these  positions. 
Suppose  three  positions  to  be  used,  and  let  the  hydrodynamic 
coefficients  be  defined  about  the  center  support  position.  The 
motion  equations  of  the  model  are 
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(72)  Iz?=  Np3  + Kr^+  Np3  + N^+  N0 

(73)  Mf  - 3)  = Yp3  + Y£$/  + Yp0  + *$f  + Y0 

where  N0  is  the  moment  about  the  central  support  position  applied 
to  the  model  through  the  shaft,  and  Y0  is  the  side  force  applied 
to  the  model  by  the  shaft.  When  the  rotation  is  about  the  center 
axis, 

(74)  H„  - -Kite  - n) 

When  the  rotation  is  about  the  forward  axis,  the  applied  moment  is 
given  by 

(75)  N0  * -Kg(f  - 72)  + YqA 

where  A is  the  distance  between  the  center  and  forward  support 
positions . 

The  angle  of  sideslip  at  the  center  support  position  is  equal 
to  yr  when  the  rotation  is  about  the  center  support.  When  the  model 
is  rotated  about  the  forward  support,  the  angle  of  sideslip  is 
given  by 

(76)  P 

Hence,  when  the  forward  support  position  is  the  center  of  rotation, 

(77)  Y0  = -mA^-  Yp(/  + Y$  - Y pf  - Y - Tj>f 

and  the  applied  moment  is 

(78)  H0  = -£ 2ty  - y2)  - A ["(mA  + Y£  + ^ 

+ (Yr  +^Yp  + Y£)?r+  Yp^] 

Equation  72  then  becomes 

(79)  [l,  - Up  - “ Bp  + A(rnA  + If  +-^-Ij)]r 

- [*r  ♦ + Tp  + ^VJr 

- [Bp  - MpV  . -KgOr-  /2) 
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A A 

P = Iz  - Np  Np  + A(nA  + Yj  + — Y$) 

-Q  = U«  + No  + — Nft  - A(Y-  + Yo  + — Yo) 


*r  "a  ‘r — "e 

p V p 


tr  + Xp  f — Xp; 


R - Np  - YpA 


Then  Eq.  79  becomes 

(81)  ?#+  o,f-  -KgO^-  y2) 

When  the  center  support  Is  the  center  of  rotation  A » 0.  The 
coefficients  for  this  case  will  be  denoted  by  Pc,  Qq,  and  Rc.  The 
coefficients  for  the  forward  support  position  will  be  denoted  by 
Pp,  Qp,  and  Rp. 

The  equation  relating)/,  5^,  and  7g  88  obtained  from  Eq.  7^ 

Is 

(82)  r - (1  + — )£ * 

Kg  Kg 

Let  the  Input  displacement  angle  7^  be 

(83)  7^  = Aj  sin  COt 
The  measured  output  Is 


/g  **  Ag  sin  (Ot  - 6) 

= Ag  cos  € sin  COt  - Ag  sin  6 cos  (Ot 


Hence,  from  Eq.  83, 


r Ki 

y = (1  + — )Ag 

*2 


*1  1 # 

cos  S A,  sin  Cot 

**»  • 


(l  + — )Ag  sin  € cos  <0t 
Kg 
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Let 

(86) 


(1  + — )A2  cos  € - — Ax  = a,  - (1  + — )Ag  sin  £ = b 


*1 


K, 


k2 


Solution  of  Eq.  8l  then  yields 

(87)  K^A^Ag  sin  s 

&)(a^  + b^) 


(88) 


R +d)2p  « K2  - K2A2 


K1  KX 

(lt^)Aa'^Al  C03ff 

a2  + b2 


Since  the  right-hand  sides  of  Eq.  87  and  88  contain  only  meas- 
urable quantities,  it  is  possible  to  determine  the  combinations  of 
hydrodynamic  coefficients  of  Eq.  80.  Thus 

(89)  Pc  = Izc  - Hp 

-Qc  » Hp  + Kg 
Rc  « h 

where  Iz?  is  the  moment  of  inertia  about  the  center  support  when 
the  model  13  rotated  about  the  center  support  position.  Rotation 
about  the  forward  support  position  yields 

(90) 

-«F  - »r  + - A(Yr  + Ip 

Rp  = Np  - A?p 

where  lzp  is  the  moment  of  inertia  about  the  center  support  posi- 
tion when  the  center  of  rotation  is  the  forward  support.  It  may 
differ  from  Izc  because  the  mass  of  the  model  may  be  changed 
slightly  in  changing  the  center  of  rotation.  The  mass  of  the 
model  when  the  forward  support  position  is  used  is  denoted  by  mp. 
With  three  support  positions  Eq.  89  and  90  yield  all  the  hydro- 
dynamic coefficients  if  it  is  assumed  that 


! 
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If  only  two  support  positions  are  used  it  is  necessary  to  estimate 
one  of  them,  and  Yg  is  the  most  convenient  to  estimate. 

Another  type  of  forced  oscillation  test  may  be  used  to  deter- 
mine the  lift  coefficient.  The  model  is  made  to  oscillate 
laterally  in  the  water  tunnel  by  means  of  a sinusoidal  force 
applied  through  a drive  spring  as  shown  in  Pig.  2. 


Let  displacement  of  the  lower  end  of  the  spring  be  £ and  the  dis- 
placement of  the  upper  end  ??.  The  force  on  the  model  is 

(91)  Y0  - " "*) 

where  K is  the  constant  of  the  spring.  The  motion  equation  of  the 
model  is 

(92)  (mV  + Yp)|3  + Ypp  + Y0  = 0 
The  angle  of  sideslip  is 

(93)  3 = -4 

V 

Hence,  using  Eq.  93 > Eq.  92  becomes 

Y*  Y 

(<*)  <*  + — )S . (-)*  + tn  - k/ 

V V 

If  the  measured  response  to  an  input, 

(95)  £ » Ax  sin  art; 
is 

(96)  c - A2  sin  (art;  - €) 


l 

i 

t 

4 

j 

I 

4 
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Eq.  $k  may  be  solved  for  Yg  and  Yg  as 


(97) 


KV  Ax 

Yq  sin  £ 

(o  A2 


KV  A2 

Yq  = — (1 cos  e)  - mV 

S 0)2  A2 


OTHER  SOURCES  OP  HYDRODYNAMIC  COEFFICIENTS 

At  the  present  time  model  tests  carried  out  in  water  tunnels 
or  towing  tanks  are  the  most  dependable  sources  of  hydrodynamic 
coefficients.  For  preliminary  design  purposes,  however,  the 
estimation  of  the  hydrodynamic  characteristics  of  a torpedo  yet  to 
be  built  is  essential.  This  problem  will  become  more  acute  when 
torpedoes  of  higher  speed  and  more  complicated  trajectories  are 
built.  It  will  then  be  necessary  to  specify  a body  shape  and  a 
tail  configuration  that  will  permit  desired  performance  while  the 
torpedo  Is  still  in  the  drawing-board  phase  of  development. 
Attempts  are  being  made  to  estimate  coefficients  on  the  basis  of 
empirical  data  and  hydrodynamic  theory.  The  aim  of  this  work  is 
to  determine  the  hydrodynamic  coefficients  given  a body  shape  and 
tail  configuration,  or,  specifying  coefficients,  to  construct  a 
body  and  tail  configuration  having  the  desired  hydrodynamic  coeffi 
cients.  Model  tests  would  then  be  only  a check  on  the  preliminary 
values  if  these  attempts  are  successful. 

The  mass  accession  terms  of  the  motion  equations  are  usually 
estimated  since  they  cannot  be  determined  from  static  or  rotating- 
arm  tests  on  models.  It  is  customary  to  evaluate  these  terms  by 
assuming  the  torpedo  to  be  an  ellipsoid  of  revolution.  The  coeffi 
cients  a^j  of  Eq.  36  are  then  given  by 

all  * klaf 

a 22  * ®33  25 

aiiJi.  3 0 

a55  - a66  a k,If 
a26  = a35  = 0 
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where  Is  the  mass  of  the  displaced  fluid  and  If  Is  the  moment 
of  Inertia  of  the  displaced  fluid  about  a minor  axis  of  the 
ellipsoid.  The  k^,  k2>  and  k'  are  Lamb’s  coefficients  (see  Ref .8). 
It  has  been  assumed  that  the  origin  is  at  the  center  of  the 
ellipsoid.  If  the  origin  is  forward  of  the  center  a distance  A 
the  coefficients  ai(j  are  given  by 


all  = klmf 
a22  ~ ^33  ~ k^f 

aH  - o 

a55  = a66  = k'lf  + k2mfA2 
a35  ~ “a26  = k2mf^ 

A check  on  the  estimated  coefficients  may  be  obtained  from 
free-fllght  tests  of  instrumented  torpedoes.  When  records  of 
pitch,  yaw,  roll,  depth,  and  control-surface  deflections  are 
obtained  from  full-scale  torpedoes  in  free  flight,  a comparison  of 
the  recorded  response  with  the  response  computed  from  the  equations 
of  motion  may  be  made.  This  avenue  of  approach  has  not  as  yet  been 
fully  exploited,  and  much  work  remains  to  be  done  before  techniques 
are  developed  for  determining  hydrodynamic  coefficients  from  free- 
fllght  records. 


FIELDS  OF  APPLICATION  OF  MOTION  EQUATIONS 


The  nonlinear  motion  equations  (Eq.  5*0  are  very  complex. 
Moreover,  in  order  that  they  be  exact,  it  is  necessary  that  the 
hydrodynamic  coefficients  of  the  left-hand  members  of  the  equations 
be  considered  as  nonlinear  functions  of  the  velocity  components. 

It  is  improbable,  therefore,  that  the  complete  equations  in  this 
form  will  ever  be  of  great  use  to  engineers . The  partially 
linearized  equations  (Eq.  55),  on  the  other  hand,  may  be  solved 
without  great  difficulty  with  the  use  of  an  analog  computer.  In 
cases  where  they  are  valid,  these  simplified  equations  can  provide 
valuable  information  about  the  trajectory  of  a torpedo.  In  most 
studies  that  have  been  made  up  to  the  present  time  the  pitch  and 
yaw  equations  (Eq.  56  and  57)  have  been  deemed  sufficient. 
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SOLUTION  OP  STEADY-STATE  EQUATIONS 

The  motion  equations  (Eq.  5k)  may  be  solved  if  a steady-state 
motion  exists,  such  as  trimmed  straight  flight  or  a steady-state 
turn.  Solution  of  the  equations  in  these  cases  is  useful  because, 
in  the  one  case,  the  trim  pitch  angle  and  elevator  setting  must  bo 
known  in  order  to  set  the  control  system  for  straight  flight;  and 
in  the  other  case,  the  change  of  depth  in  a steady-state  turn  may 
be  determined. 

Trim  Flight 

If  the  torpedo  is  moving  in  straight  trimmed  flight 
u=^=]/-=^=^=©  = 0 = p = p = a = d°I>  = O 

and 

8 = a 

The  motion  equations  (Eq.  5k)  become 

(98)  XQ  + Tx  - (W  - B)a  + ^u  = 0 

Zaa  + V/e  + O'  - »)  - 0 


Maa  + Ve^e  + BxB  " WzGa  “ 0 


Solving  for  the  trim  values  of  u,  a,  and  gives 

-ZaBXl  • (W  - B)(Ma  - Wz5) 

n *yve  - veov  - 

-Z^  Bxb  + McT  (W  - B) 
a*  S e 

V,(«a  - WzG>  - «(feV 
X0  + (W  - B)ttl 

Ml  = + 

-Xu  X* 
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Steady-State  Trim 

If  the  rudder  is  given  a steady  deflection  when  the  torpedo 
is  in  straight  trimmed  flight,  a rolling  moment  will  be  produced 
because  of  the  centrifugal  force  acting  on  the  center  of  gravity 
of  the  torpedo,  which  lies  below  the  longitudinal  axis.  Because 
of  the  resulting  roll  the  rudders  will  cause  the  torpedo  to  spiral 
downward.  However,  the  change  in  depth  will  cause  the  control 
system  to  function,  giving  an  up  elevator  which  will  tend  to  re- 
duce the  depth  error.  When  a steady  state  is  attained,  the 
torpedo  will  circle  at  a constant  depth  error  with  constant  angles 
of  roll,  pitch,  sideslip  and  attack,  and.  a constant  elevator 
deflection  different  from  the  trim  value.  Under  these  conditions 


Y = constant 

and 

• • «••#,  «•  • • 
u = © = (p=Q=j/=.(p  = aJ  = (3  = 0 

The  motion  equations  (Eq.  5*0  become 

(99)  X0  + Xuu  + Tx  - (W  - B)0  = miV(a  sin  <p  + p cos  tp)y 

+ y2(-Z<j  sin2  <p  + Yr  cos2  <p) 
- mz oj£&  cos  (p 

Ypp  - Yp y&+  Yp^cos  p + Y^c/’p  + (W  - B)  sin  <p 

= mLU y cos  p + mzgjK2  cos  p sin  (p  + m£Ja.yr& 

- Z^9  sin  <p 

Z^a  + ZqY  sin  <p  + Z<f  <?e  + (W  - B)  cos 

= m^V^fop  - mzjj^2®2  - mzg y2  sin20>  + Yj^2®  cos  (p 

- mjU^  sin  p 
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(99)  Kp0  - Kpjfe  + Y^jr  cos  <p  + Vp</r  - Wzg,  sin  <p 
Contd . 

= (Jz  - Jy)^  cos  <p  sin  <p  + Ixz^%  sin  <p 
+ VjKY^  + Zq)(p  sin  <p  + a cos  <p ) - mzGU^  cos  <p 
- mzQVa^ 

Maa  + Mq$y  sin  (p  + ^ + Bxg  cos  (p  - WzG# 

- - J*)r^  COS  {P  + i^2  - V2  cos2  p 

- Nppjfe  + mzGV^(a  sin  q>  + 0 cos  <p) 


Np3  - Np y9  + N vf  cos  <p  + N <f  <fv  - Bxg  sin  <p 

- (Jx  - Jy )p20  sin  <p  + M^a fs  + Ixz^  cos  <p  sin  <p 
0 - a cos  <p  - 3 sin  <p 

These  equations  are  to  be  solved  for  y,  (p,  0,  a,  9,  u,  and  <f&.  A 
numerical  method  of  successive  approximation  is  probably  most  con- 
venient to  use.  A first  approximation  may  be  obtained  by  line- 
arizing Eq.  99  in  all  variables  except  jp.  A first  approximation 
to  <p  is  given  by 

_ va / 

(100)  (p  = tan"1  — . 

g 

The  following  equations  then  yield  first  approximation  to  y and  0. 

(101)  Yp0  + Yry  cos  <p  + X + (w  “ B)  sin  <f  = m^V^  cos  <p 

Np0  + N py  cos  <p  + N</»  d*T  - Bxb  sin  (p  = 0 


Using  the  first  approximations  to  yr  and  0 obtained  from  Eq.  101 
Eq.  102  may  be  used  to  yield  first  approximations  to  a.,  and 


3k 


I 

1 


t 


$ 

: 


4 


I 


i 


i 


(102)  Zaa  + Zqft  sin  <p  + Z^^e  + (W  - B)  cos  <f  = - m -gif  sin  (p 
Map,  + Mq^  sin  (p  +■  M <Pe^e  + Bxg  tcos  (p  - Wzq©  = 0 
0 = a cos  (p  - p sin  (p 

If  more  accurate  solutions  are  desired  a numerical  method  of  suc- 
cessive approximations  may  be  used  starting  with  the  first 
approximations  given  by  Eq.  100,  101,  and  102  (see,  for  example, 
Numerical  Calculus,  W.  E.  Milne,  Princeton  University  Press, 1949 ) . 

Depth  error  in  the  turn  may  be  calculated  by  considering  the 
depth  error  necessary  to  yield  the  elevator  setting  given  by 
Eq.  99- 


STABILITY  OP  CONTROLLED  TORPEDOES 

The  behavior  of  a torpedo  in  the  water  is  a function  not  only 
of  its  hydrodynamic  characteristics  but  also  of  its  internal 
control  system.  The  complete  system  must  be  considered  before  it 
can  be  decided  whether  a torpedo  is  capable  of  the  performance 
that  is  required . Study  of  the  complete  system,  comprised  of 
hydrodynamic  characteristics  and  control  system,  is  usually  termed 
"stability  analysis" . It  Is  not  the  purpose  of  this  report  to 
discuss  all  the  methods  by  which  stability  analyses  may  be  under- 
taken, but  it  seems  appropriate  to  describe  the  manner  in  which 
the  motion  equations  enter  into  the  problem.  A control  system  for 
a torpedo  contains  devices  that  can  detect  the  position  or  attitude 
of  the  torpedo  or  their  rates  of  change.  Signals  from  these  de- 
vices are  used  to  control  the  action  of  elevators  or  rudders  which 
produce  changes  In  the  trajectory  of  the  torpedo.  Thus  the  con- 
trolled torpedo  constitutes  a feed-back  system,  or  servomechanism, 
and  standard  analysis  techniques  from  the  theory  of  servomecha- 
nisms are  applicable  (see  Ref.  9)* 

A simple  example  will  be  used  to  show  how  a control  system 
may  be  analyzed.  Assume  that  the  trajectory  of  a torpedo  is  in  a 
horizontal  plane  and  that  the  torpedo  does  not  roll.  The  motion 
equations  are  given  (Eq.  56)  as 

(103)  - mTVp  » YgfJ  + (Yr  - mLV)^  + Yr^  + y/v 

(104)  Jzf  = Npp  + + Np3  + ^4 


NAVOBD  REPORT  2090 


The  rudder  Is  controlled  by  a signal  obtained  from  a device  sen- 
sitive to  direction,  such  as  a gyroscope.  It  will  be  assumed  that 
the  rudder  deflection  is  proportional  to  the  difference  between 
the  torpedo  direction  yr  and  a reference  direction  • A time  lag 
Tg  is  Introduced  in  the  equation  to  represent  the  delay  in  the 
actuator  that  operates  the  rudder. 

Thus 

d 

(105)  (1  + Tg )cTt  = Kj irty  - ys) 

dt 

Equations  103,  104,  and  105  may  be  solved  for  y in  terms  of 
V's  • Usually  y/s  Is  a constant  direction.  Depending  on  Ky  and  the 
hydrodynamics,  the  motion  of  the  torpedo  will  be  either  stable  or 
unstable.  It  is  said  to  be  stable  If  transient  oscillations  are 
eventually  damped  out  and  the  torpedo  assumes  the  direction  ys  • 

The  motion  is  unstable  if  the  oscillations  continually  Increase  In 
amplitude . It  is  the  aim  of  the  designer  to  choose  a YLy  that  will 
result  in  a fast  well-damped  response  to  a disturbance. 

Even  a simple  system  such  as  that  of  the  example  cited  above 
requires  long  tedious  computation  if  it  is  to  be  solved  without 
recourse  to  mechanical  or  electrical  computers.  A great  saving  In 
time  and  effort  is  achieved  by  the  use  of  such  aids  to  computation. 
This  type  of  problem  is  particularly  amenable  to  solution  with  the 
use  of  an  analog  computer  such  as  the  REAC  (Ref.  10).  Appendix  6 
shows  how  Eq.  103,  104,  and  105  are  solved  on  the  REAC  at  the 
U.  S.  Naval  Ordnance  Test  Station.  The  REAC  is  capable  of  solving 
much  more  complicated  problems,  but  the  principle  of  operation  is 
demonstrated  rather  well  by  this  example. 
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aU 

B 

F 

2l 

E2 

f 

g 

Gjj 

Of 

G 

Sb 

Sf 

H 

Hx,  Hy,  Hz 
*x>  *y>  IZ 
Iyz>  Ixz>  Ixy 
Jx>  ^y>  Jz 

k 


Appendix  A 
NOMENCLATURE 

Components  of  apparent  mass  tensor 
Buoyancy  of  torpedo 
Force  acting  on  torpedo 

Force  on  torpedo  predicted  from  potential  flow 

l - 1 1 

Force  per  unit  volume  acting  on  fluid 
Acceleration  of  gravity 
Linear  momentum  of  torpedo 
Linear  momentum  of  fluid 

Sf  + Sb 

Angular  momentum  of  torpedo 
Angular  momentum  of  fluid 

Sf  + Sb 

Components  of  H in  body  coordinates 

Moments  of  inertia  about  x,  y,  z axes,  respectively 

Products  of  inertia 

Apparent  moments  of  inertia  about  x,  y,  z axes, 
respectively 

Viscosity  coefficient  of  fluid 
Components  of  Lg  in  body  coordinates 
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j 

I 

i 

i 

j 

j' 


L Moment  acting  on  torpedo 
Li  Moment  on  torpedo  predicted  from  potential  flow 
Lg  L - Lx 
m Mass  of  torpedo 

mi,  Apparent  longitudinal  mass  of  torpedo 
mip  Apparent  transverse  mass  of  torpedo 
n Unit  normal  to  surface 


p,  q,  r Components  of  uj  in  body  coordinates 
P Pressure  acting  on  fluid 

Velocity  field 
T Thrust  of  torpedo 

Tx  Magnitude  of  thrust  of  torpedo 
Tjj  Kinetic  energy  of  torpedo  body 

Tf  Kinetic  energy  of  fluid 

T Tb  + Tf 

U,  v,  w Components  of  V in  body  coordinates 
V Velocity  of  torpedo 

W Weight  of  torpedo 

x , y,  2 Body-fixed  coordinates 

x0,  yo>  zo  Space-fixed  coordinates 

xB,  yB,  2B  Components  of  buoyancy  moment  arm 
Xq,  y q,  zq  Body  coordinates  of  torpedo  c.  g. 

X2,  z2  Components  of  Pg  in  body  coordinates 
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W, 


V,e, 


a Angle  of  attack  of  torpedo 

0 Angle  of  sideslip  of  torpedo 

J'q  Elevator  deflection 

<fr  Rudder  deflection 

V Direction  cosines  of  n 

Density  of  fluid 
0^  Velocity  potential  function 

(p  Inertial  reference  angles  of  torpedo 

*0  Angular  velocity  of  torpedo 
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Appendix  B 

PORCE  AND  MOMENT  ON  TORPEDO  PREDICTED  PROM  IDEAL  FLUID 


Consider  the  torpedo  to  be  immersed  in  an  ideal  fluid  of 
infinite  extent,  the  fluid  being  at  rest  at  infinity.  Let  the 
velocity  field  be  denoted  by  The  total  linear  momentum  of  the 
fluid  is  given  by 

(106)  if  =fTPl  <iT 

where  the  integral  is  over  the  entire  volume  of  the  fluid.  The 
momentum  of  the  fluid  must  have  a finite  value  since  otherwise  it 
would  be  Implied  that  an  Infinite  momentum  had  been  Imparted  to 
the  fluid  by  finite  forces  exerted  for  a finite  time  by  the  tor- 
pedo, and  this  is  Impossible.  Now  consider  an  element  of  the 
fluid  occupying  the  volume  d T.  By  Newton's1  second  law  the  force 
dP  acting  on  the  element  Is 


d 

(107)  dP  = — (/3a  dT) 

dt 

The  total  force  acting  on  the  fluid  is  obtained  by  integrating 
over  the  whole  fluid.  Then,  if  Ff  denote  the  total  force  acting 
on  the  fluid. 


d 

(108)  Pf  */  — (p£  d7) 

T dt 

d 

■ — Lpsl  iT 

dt  " 
dOf 

SB  " 

dt 


The  total  force  on  the  fluid  is  that  exerted  by  the  torpedo.  Hence, 
by  Newton's  third  law  the  force  Fj.  on  the  torpedo  is 


h 9 “*f 


dt 


(109) 
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In  a similar  manner  the  net  moment  acting  on  the  fluid  about 
a point  in  Inertial  space  may  be  shown  to  be  equal  to  the  inertial 
time  rate  of  change  of  the  total  angular  momentum  of  the  fluid. 
Then,  also, 

dHf 

(110)  

dt 

where 

(111)  Hf  =f^p  s ar 

S being  the  radius  vector  of  the  fluid  element  /OdT. 
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Appendix  C 

FORCE  AMD  MOMENT  REFERRED  TO  A MOVING  COORDINATE  SYSTEM 


The  laws  of  dynamics  state  that  the  external  force  acting  on 
a system  Is  equal  to  the  time  rate  of  change  of  linear  momentum 
of  the  system,  and  that  the  moment  of  this  force  about  a point 
fixed  In  space  Is  equal  to  the  time  rate  of  change  of  the  angular 
momentum  about  this  point.  Let  Oxyz  be  the  body-fixed  coordinate 
system  of  the  torpedo,  and  let  the  velocity  V of  0 have  the  com- 
ponents U,  v,  w on  these  axes . The  angular  velocity  U)  about  0 
has  the  components  p,  q,  r.  Let  G be  the  linear  momentum  of  the 
system  and  let  H1  be  the  angular  momentum  about  a point  O'  fixed 
In  Inertial  space.  Let  H be  the  angular  momentum  about  0.  The 
time  rate  of  change  of  a quantity  seen  from  O'  will  be  denoted 
by  the  operator  d/dt.  The  time  rate  of  change  as  seen  from  the 
moving  coordinate  system  Oxyz  will  be  denoted  by  a dot  placed  over 
a symbol. 

The  force  acting  on  the  system  is  given  by  (see  Ref.  6) 

dG 

(112)  F = — . G + a)  X G 

dt 

The  angular  momentum  H1  about  0>  and  the  angular  momentum  H 
about  0 are  related  by 

(113)  H1  - H + S x G 

where  S is  the  radius -vector  from  0*  to  0.  Thus  the  moment  L1 
about  O'  is  given  by 

dH'  dH 

(114)  L«  = — • = — + SXF  + VXG 

dt  dt 

since  V «*  dS/dt.  The  moment  L about  0 Is  given  by 
(H5)  L = L«  - S X F 

Hence, 

(116)  L = H+  dJXH  + VXG 
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Appendix  D 

DERIVATION  OP  MOMENTUM  PROM  KINETIC  ENERGY 


The  kinetic  energy  of  an  ideal  fluid  is  given  in  Eq.  31  as 

1 r dtp 

(117)  T = —pl0 — dS 

2 X dn 

where  the  integral  is  taken  over  the  surface  of  the  torpedo.  The 
potential  function  has  the  *'orm 

(118)  <p  = + vfy  + w + p^  + q (fo  + v<j>6 

The  velocity  field  c[  of  the  fluid  is  obtained  from  the  velocity 
potential  as 

(119)  a ■ ’ V0> 

The  total  linear  momentum  of  the  fluid  is  obtained  by  integrating 
the  momenta  of  the  mass  elements.  Thus 

(120)  Gf  = 

By  means  of  the  divergence  theorem  the  volume  integration  may  be 
expressed  as  an  integration  over  the  torpedo  surface  as 

(121)  £f  ^ J1  dS 

where  n is  defined  as  in  the  section  "Kinetic  Energy  of  Ideal 
Fluid.'”’  The  component  of  G in  the  x direction  (using  Eq.  2k)  is 

(122)  Gf 

x 

in  view  of  Eq.  25-  Using  Eq. 

(123)  Gf^  = al]G  + al2v 


•pji?*  ® 


118,  Eq.  122  becomes 
+ a13w  + a^p  + a15q  + a^ 
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It  Is  seen,  then,  that 

&Tf 

(124)  Gf  = 

* &U 

It  may  be  shown  In  a similar  manner  that  the  other  components 
of  Gf  are 

wf 

(125)  Gf  = 


fc*f 

Gf  = - — 

2 *W 

The  total  angular  momentum  of  the  fluid  is  given  by 


(126) 


5f 


■ ~P  Is  £ 


d T 


where  S is  the  radius  vector  of  the  fluid  element  p dT.  Now 
(12?)  V<p  x S « V X pS  -g)\j  x S 

=v  x <ps 


since  7 x s a 0 

By  the  divergence  theorem,  then,  the  volume  Integral  of 
Eq.  126  may  be  transformed  to  an  integral  over  the  surface  of  the 
torpedo  as 

(128)  Hf  = -pfs<pn  X S dS 

The  component  of  Hf  in  the  x direction  is  therefore 

(129)  Hf  = -p[  (PiJJz  - Vf)  dS 

x JS 


with  reference  to  Eq.  27.  Using  Eq.  118  and  32,  Eq.  129  may  be 
written 


a^U  + a2j+v  + aj^w  + a^yp  + aj^q  + a^gr 
dTf 


(130) 
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Similarly, 

(131)  Hf 

l7 


The  components  of  momentum  of  the  torpedo  body  may  be  deter- 
mined from  the  kinetic  energy  of  the  body.  The  linear  momenturf 
is  defined  by 

(132)  -[  mi(V  + to  x r±) 
and  the  angular  momentum  by 

(133)  Hb  = •£  m^  x (V  + 0)  x r±) 

where  the  summation  is  over  all  particles  m.^  of  the  torpedo,  and 
rj  is  the  radius  vector  from  the  origin  to  n^. 

Since 

1 

(134)  T b = - T JHiCV  + CO  X 

2 

iTb  _ ZV 

(135)  ■ ) m±(V  + CO  x r^X — + — x r*  + W x ) 

6U  dU  bV  " 6U 

= l mi(V  + 03  X Zi)-1 
= Gbx 

6Tb  iTb 

Similarly,  Gv  and  — = Ck 

bv  y z 

Differentiating  Tb  with  respect  to  p yields 

6 Tb  6V  buj  6r^ 

(136)  . £ m^V  + OJ  X rjX~  + — X + {0  X ) 

6p  " ~ " 3p  bp  ” “ip 


6Tf 

6q 

6Tf 

6r 
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(136) 
Contd . 


- Y ®i(V  + w x riXi  x r±) 
~ Y * X (V  + (O  x r^)*i 


Similarly, 


dTb  iTb 

S~vir-^ 


Consequently,  since  T = + Tf,  G = Gf  + G^,  and  H = Hf, 


(137) 


dT  dT 

G as 1 + l +—  ic 

dU  dv  dw 

dT  dT  dT 

H i + J + k 

dp  dq  dr 
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Appendix  E 

RESOLUTION  OP  BUOYANCY  AND  GRAVITY  FORCES  AND 
MOMENTS  ONTO  BODY  COORDINATES 


The  orientation  of  a torpedo  with  respect  to  a fixed  inertial 
frame  of  reference  is  specified  by  an  angle  of  yaw  y,  an  angle  of 
pitch  O,  and  an  angle  of  roll  cp.  Let  the  inertial  axes  be  x0, 

Yo>  zo>  and  let  the  unit  vectors  in  the  direction  of  these  axes  be 
io>  respectively.  The  body  axes  at  the  start  coincide 
with  the  inertial  axes  and  are  then  rotated  about  the  zQ  axis 
through  an  angle  ^to  coincide  with  axes  xj_,  y^,  zj.  Let  the  unit 
vectors  in  the  directions  of  these  axes  be  3^,  ^l,  k^.  Then 

(138)  ii e io cos  y + io sin  y 

ll  - -io  sin  y + lo  cos  y 


Now  the  body  axes  are  rotated  about  the  y^  axis  through  an 
angle  & to  coincide  with  axes  x2,  Y&  z2'  ^tting  the  unit 
vectors  in  the  directions  of  these  axes  be  ig,  jg,  kg, 

(139)  ig  = i^  cos  6 - kj  sin  6 


^•2  “ ^1 


Finally,  the  body  axes  are  rotated  around  the  x2  axis  through 
an  angle  <p.  Letting  the  body  axes  be  x,  y,  z and  the  unit  vectors 
in  the  directions  of  these  axes  be  i,  k,  respectively,  yields 

(140)  i = ig 

= Jg  cos  <p+  kg  sin  <p 
k = -Jg  sin  (p+  kg  cos  (p 
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Combining  these  transformations  gives 

(l4l)  i = Iq  cos  Y cos  & + Ao  3ln  Y cos  ® “ £o  sln  © 

1 **  lo(cos  Y sin  0 sin  (p  - sin  ^cos  <p) 

+ 2o(sin  Y sin  & sin  f + cos  Y cos  fft 
+ kg  cos  & sin  (p 

k a ig(sin  y sin  (f  + cos  Y sin  £ 003  p) 

+ ^(si nY sin  £ cos  ¥ ~ cos  Y sln  $ 

+ ^ cos  9 cos  p 
or 

(lk2)  Iq  * i cos  y cos  9 + j(cos  y sin  9 sin  <jP  - sin  ^cos^>) 

+ k(sln  Y sin  ¥ + 003  Y 3in  ^ cos  $ 

Jo  a i sin  Y cos  & + j(cos  Y 003  ¥ + sin  Y sin  & sln  0 

+ k(sin  Y sin  & 003  ¥ " cos  Y 3in 

ko  » -i  sin  9 + J cos  & sin  (/>  + k cos  6 cos  cf> 

The  weight  and  buoyancy  forces  act  in  the  z0  direction. 

Hence  the  vector  representing  the  gravity  buoyancy  force  is 

A 

(l43)  (W  - B)ko  a (w  - B)(-i  sin  & + J sin  ^ cos  9+  k cos  ^cos  9) 

The  components  of  the  gravity-buoyancy  force  may  be  read 
from  this  expression. 

If  £g  is  the  radius  vector  from  the  origin  of  the  body  coor- 
dinates to  the  center  of  buoyancy,  the  moment  about  the  origin 
caused  by  the  buoyancy  force  is 

W)  £g  x (-Bkg) 

Since  r a ix_  the  moment  due  to  buoyancy  is 
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(1^5)  JBxg  cos  cp  co s 9 - kB*B  sin  cp  cos  & 

Letting  rg,  be  the  radius  vector  from  the  origin  to  the  center 
of  gravity  of  the  torpedo,  the  moment  about  the  origin  caused  by 
gravity  is 

(146)  £q  x 

Since 


£G  = ZG* 


the  gravity  moment  is 

(147)  - iWz(j  sin  cp  cos  6 - jWzg  sin  9 


k9 
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Appendix  F 

RESOLUTION  OF  ANGULAR  VELOCITY  TO  INERTIAL  COORDINATES 


The  angular  velocity  of  the  torpedo  has  the  components  p,  q, 
r In  body  coordinates.  It  is  required  that  the  angular  velocity 


be  expressed  in  Inertial  coordinates  0,0).  ^ ls  noted  that 

the  z0  axis,  6 . 


^"is  the  angular  velocity  about  the  z0  axis, 9 is  the  angular 
velocity  about  the  y^  axis,  and  (f  is  the  angular  velocity  about  the 
x axis.  Hence 


(148)  <0  = ip  + Jq+lcr  = i^  + + kjj/ 

Using  the  relations  between  the  unit  vectors  of  Appendix  E, 

p s (f>  - y sin  9 

• 0 

q = <?  cos  (f  + y sin  <5^  cos  0 

• • 

r = -0sin  (p  +Y  cos  ¥ cos  & 


(149) 
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Appendix  0 


REAC  ANALOO  OP  TORPEDO  SYSTEM 


The  system  to  be  solved  on  the  REAC  is  described  by  the  sys- 
tem of  equations 


(150) 


- mTVDp  = Ypp  + (Yr  - mLV)j/+  Ypf  + Y j 


3ZW=  NgP  + NrJ/+  NgDp  + ^/T 
(1  + TfD)cfr  = s) 

where  D represents  the  operation  of  differentiation  with  respect 
to  time  t.  Equations  150  are  more  readily  analoged  if  they  are 
transformed  as  follows: 

(151)  -(mijiV  + Yj  ~ )Dp  * (Yjj  + Yp  )p  + (Yr  - mLV  + Y£  — 


N<fr 

+ (Y<fr  + Y* )(TP 

J, 


Yr  - mLV  ^ 
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The  REAC  analog  of  this  system  is  shown  in  Pig.  3*  The 
symbols  used  have  the  following  meaning: 


= Summing  Amplifier 


= Integrating  Amplifier 


= Gain  of  amplifier  number  i 
with  input  number  J 


= Voltage  Divider 
(Potentiometer) 

= Setting  of  potentiometer 
number  i 


The  amplifier  inputs  have  available  gains  of  1,  and  10. 
The  output  is  inverted  in  sign. 

The  equations  of  the  analog  are 


(152)  ex  e2  e3 

el  ■ a1.3k1  ~ + a1.2a6k2  ZT  " a1.1k3  “ 

D D EF 

e2  el  e3 

e2  = “a2.3k4  “ + a2.1a5k5  — ■ a2.2k6 


e3  ei  e4 

e3  = -a3.3k7  y + a3.1a4k8 ~2~  A3-2ZT 
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where  D is  the  operator  representing  differentiation  with  respect 
to  the  computer  time  base  t.  Let  t = nt.  Then  D = nD.  The 
voltages  e^,  e2>  e^  are  the  electrical  analogs  of  the  variables  y, 
P,  <fr,  respectively,  of  the  physical  system.  Let 

(153)  ~y 

a2e2  * 0 
a3e3  ~ 

akeb  = fa 

Equations  152  may  then  be  written 

(154)  tjr.  -nA  K/  + nA12A6K2|3  - — n ^ 

a2  a3 

a2  . a2 

Dp  = -nA0  _K4P  + — nAg  nAg  2Kg<£ 

al  a3 

D<r  - -"A3.3KTt^  +—  »S.iV6f-— »*3.aV. 

r aj  a^ 

By  comparing  Eq.  l$k  with  Eq.  151  it  is  seen  that  the  two 
systems  are  equivalent  if 


a. 


NP  " Np 


*A1. 2*6*2 


mipV 


a + 


Jz  + Np 


Hr  - Ip' 


y. 

. xr 


nA1.3Kl 


m^V 

?r  - mLV 

mijV 


Jz  ♦ 


xr 

m^V 


9* 
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— nAl . 1k3  ■ 
a3 


K <f  - 

r 


UlrpV 


Jz 


nA2.3Kk  ~ 


a2 

— nA2.1A5K5  * 

al 


Hr. 


y3  + Yr‘ 


fflipV  + Yp 


Np 


•Yr  - mLV  + Y*. 


m^V  + 


N(3 


Nr 

JZ 


a2 

— nA2.2^6 
a3 


Vr  + Yf 


nupV  + Y^ 


1 

nA3-3K7  =“ 
T1 


a3  p 

— n 2a3. 1*4*8  = — 
al  Tl 


— nA3<2 
a4 
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Amplifier  gains,  potentiometer  settings,  time  base  change  n,  and 
the  ratios  a^/^,  a^/a^,  ai/a^,  are  chosen  in  such  a manner  that 
the  above  equations  are  satisfied.  The  response  to  an  initial  yaw 
error,  obtained  by  placing  the  appropriate  initial  condition  on 
amplifier  A4,  may  be  found. 
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